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ABSTRACT 


Design  considerations  for  a  heterodyne  spatial  tracking  system  utilizing  pupil-plane  processing 
techniques  and  its  advantages  over  traditional  focal-plane  processing  are  described.  Noise  performance 
bounds,  optimal  and  suboptimal  local  oscillator  distributions,  pull-in  performance,  and  applications  other 
than  spatial  tracking  are  discussed.  Experimental  verification  of  a  one-axis  closed-loop  tracking  system 
is  presented. 
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OPTICAL  SPATIAL  TRACKING  USING  COHERENT 
DETECTION  IN  THE  PUPIL  PLANE 


1.  INTRODUCTION 

Almost  all  radar  and  communications  systems  require  stabilization  to  a  fraction  of  an  antenna 
beamwidth.  At  optical  frequencies,  where  antenna  beamwidths  can  be  extremely  narrow  (~  pradians), 
spatial  tracking  issues  play  a  dominant  role  in  system  design  and  performance.  Spatial  tracking  systems 
based  on  coherent  (heterodyne,  homodyne)  or  noncoherent  (direct)  detection  can  be  used.  Although 
direct  detection  systems  are  the  most  commonly  used  and  analyzed,  coherent  detection  is  the  preferred 
choice  in  many  instances.  Such  instances  arise  when  increased  sensitivity  is  needed,  especially  in  the 
presence  of  background,  or  in  coherent  radar  and  communication  systems  where  it  is  natural  to  integrate 
a  coherent  tracker  and  a  coherent  receiver. 

Spatial  tracking  systems  which  employ  coherent  detection  have  been  the  subject  of  many  stud¬ 
ies  [1-6].  These  studies  have  focused  on  optical  detection  in  the  focal-plane,  i.e.,  focal-plane  processing. 
In  the  focal  plane  there  is  a  linear  relationship  between  focal  spot  lateral  displacement  and  incidence 
angle  (see  Figure  1).  A  tracking  error  sensor  can  be  formed  by  placing  a  quadrant  detector  in  the  focal 
plane  and  processing  the  outputs  to  extract  the  tracking  information  that  is  contained  in  the  relative 
amplitudes  of  the  four  detector  outputs.  Alternatively  optical  detection  may  be  performed  in  the  pupil 
plane  [see  Figures  2(a)  and  2(b)].  In  this  case  a  tracking  sensor  can  be  formed  by  placing  a  quadrant 
detector  in  the  pupil  plane  and  processing  the  outputs.  However,  in  pupil-plane  processing  the  tracking 
information  is  contained  in  the  relative  phases  of  the  four  detector  outputs.  Note  that  although  focal- 
plane  processing  techniques  are  applicable  to  both  heterodyne  and  direct  detection,  pupil-plane  process¬ 
ing  is  only  possible  with  coherent  detection  schemes. 
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Figure  2(b).  Pupil-plane  processing  (amplitude  processing). 
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It  is  possible  to  show  that  for  arbitrary  detector  arrays  and  electrical  processing  there  is  no  theoreti¬ 
cal  advantage,  in  the  linear  region  of  the  angle  discriminator  or  otherwise,  by  either  processing  tech¬ 
nique.  That  is  because  coherent  detection  and  all  passive  optical  systems  are  invertible  linear  transfor¬ 
mations.  Equivalently,  there  is  no  information  lost  as  the  signal  field  propagates  from  the  pupil  plane  to 
the  focal  plane.  However,  when  implementation  issues  (such  as  receiver  complexity)  are  considered,  a 
number  of  attractive  features  of  pupil-plane  processing  become  apparent.  For  instance,  in  order  to 
overcome  the  power  loss  and  discriminator  distortions  that  result  from  the  dead-zone  (between  1  to  5  pm 
in  size  for  a  good  device)  associated  with  focal-plane  quadrant  detectors  or  image-splitters  one  must  use 
long  focal  length  lenses  (large  f  numbers)  or  complex  telephoto  lenses  (see  Figure  3).  The  resultant 
system  is  usually  large  and  alignment-sensitive.  Pupil-plane  systems  are  far  less  sensitive  to  this  loss 
since  in  most  cases  the  pupil  size  is  naturally  larger  than  the  dead  zone  and  since  pupil-plane  dead  zones 
will  not  distort  the  discriminator.  As  another  example,  high-bandwidth  quadrant  detectors  with  narrow 
dead  zone  and  low  cross  talk  for  focal-plane  systems  are  difficult  to  obtain,  which  usually  results  in 
systems  which  employ  large,  bulky,  alignment-sensitive  image-splitting  techniques.  Although  these  same 
techniques  can  also  be  used  in  the  pupil  plane,  an  alternative  is  to  use  four  discrete  detectors  and  a  four- 
element  lens  to  obtain  a  compact  high  performance  detector.  As  seen  in  Figure  4,  the  dead  zone  and 
cross  talk  between  detectors  is  not  a  concern,  and  the  dead  zone  within  the  focusing  lens  can  easily  be 
made  small  compared  to  the  signal  field  diameter. 

The  following  sections  will  discuss  pupil-plane  processing  issues  in  more  depth.  The  next  section 
will  present  unbiased  Cram^r-Rao  bounds  on  the  noise  performance  of  a  pupil-plane  tracking  system  and 
compare  them  with  the  focal  plane.  Section  3  will  investigate  suboptimal  local  oscillator  distributions 
and  compare  their  performance  with  the  optimal  distribution.  In  Section  4,  a  comparison  of  pull-in 
performance  of  focal  plane  and  pupil-plane  tracking  systems  will  be  inferred  from  their  off-axis  discrimi¬ 
nator  characteristics.  Section  5  will  briefly  investigate  some  potential  applications  of  pupil-plane  pro¬ 
cessing  using  higher  order  arrays.  Experimental  results  for  a  one-axis  closed-loop  tracking  system  are 
presented  in  Section  6. 
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Figure  3.  Dead  zone  consideration 
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Figure  4.  Pupil-plane  tracking  sensor 
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2.  PERFORMANCE  BOUNDS 


In  this  section  we  develop  bounds  on  the  noise  equivalent  angle  (NEA)  of  a  heterodyne  spatial 
tracking  system  operating  in  the  pupil  plane.  A  pupil  plane  can  be  considered  a  perfect  image  of  the 
entrance  aperture  of  an  optic;'  .descope,  or  alternatively,  a  plane  in  which  the  power  distribution  is 
constant  with  angle.  Let  the  signal  and  local  oscillator  fields  in  the  pupil  plane  be  represented  by 


S  =  JP~s  S(x,y)  e 


-  j  k  (  X0x 


+  v6y) 


j  co  t 


(1) 


L  =  L(x,y)  e'  J  ‘  ®  +  461 1 1  <2> 

where  S  is  the  signal  field  distribution,  L  is  the  local  oscillator  field  distribution,  Ps  is  the  power  in  the 
signal  field,  Pj  is  the  power  in  the  local  oscillator  field,  k  =  (2tt/X)  is  the  wave  number,  K  is  the  wave¬ 
length  b>  =  (2irf)  is  the  optical  frequency,  Aro  is  the  I F  frequency,  t  is  time,  x  and  y  are  the  two  spatial 
parameters,  and  6  and  0  are  the  two  angular  tracking  errors  (assumed  small).  S(x,y)  and  L(x,y)  are  the 
complex  amplitudes  of  the  signal  and  local  oscillator  fields,  respectively.  S  and  L  are  assumed  to  be 
nominally  propagating  in  the  z  direction.  They  have  been  normalized  so  that 

oo  oo 

1  =  J I  S(x,y)  I2  dx  dy  =  Jl  L(x,y)  I2  dx  dy  (3) 

-OO  -oo 

For  simplicity  we  assume  that  S(x,y)  and  L(x,y)  are  symmetric  about  the  x  and  y  axes.  Further,  we  will 
assume  that  the  detector  array  is  an  infinite  area  quadrant  detector.  Note  that  the  performance  bounds  are 
valid  for  a  detector  array  with  an  infinite  number  of  detector  elements  since  L(x,y)  can  be  arbitrarily 
shaped  and  we  assume  linear  processing.  That  is,  a  given  method  of  linear  post  detection  processing 
(multiplying,  adding,  and  subtracting)  on  an  infinite  detector  array  can  also  be  done  with  an  infinite  area 
quadrant  detector  and  the  proper  shaping  of  the  local  oscillator  field  distribution  (magnitude  and  phase). 

If  the  sum  is  generated  by  the  sum  of  all  four  quadrants  and  the  azimuth  angular  error  signal  is 
generated  by  the  sum  of  the  right  two  quadrants  minus  the  sum  of  the  left  two  quadrants,  then 


Sum  =  /2PS  m  (0x,0y)  cos(Acot)  +  wSum(  t ) 

(4) 

Az  =  72PS  )  0X  sin(Atot)  +  wAz(  t ) 

(5) 
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where 


m(ex.9y)  =  — 


OOOO 

||  S(x,y)  L(x.y)  cos(kx0x)  cos(ky0y)  dx  dy 


|| !  L(x,y)  r  dx  dy 


oo 


Kaz^’0>) 


CWOO 

2  ||  S(x,y)  L(x,y)  sinc(kx9x/Ji)  cos(ky0y)  dx  dy 

00 _ _ 


JJl  L(x,y)  I2  dx  dy 


oo 


(6) 


(7) 


and  wSum(  t )  and  wAz(  t )  are  the  spectrally  white  Gaussian  noise  in  the  sum  and  difference  channel, 
respectively.  The  double-sided  spectral  height  of  the  noise  is  equal  to  \\v/(2t\),  where  h  is  Planck’s 
constant,  v  is  the  optical  frequency,  and  t]  is  the  quantum  efficiency.  The  magnitude  squared  of  m(0)  is 
often  referred  to  as  the  sum  channel  mode-matching  efficiency.  By  analogy,  the  magnitude  squared  of 
KAz(0)  could  be  referred  to  as  the  azimuth  channel  mode  matching  efficiency.  However,  it  is  more 
commonly  known  as  the  azimuth  channel  discriminator  gain. 

By  taking  advantage  of  the  symmetry  of  the  sum  and  difference  channels,  an  equation  describing 
the  normalized  signal  out  of  the  right  two  quadrants  can  be  derived. 1  This  equation  after  normalization  to 
the  same  noise  spectral  density  as  above,  which  illustrates  the  dependence  of  the  phase  of  the  I F  signal 
on  tracking  error,  is  given  by 


AZ+  =  y2iy 


m(0x,0v)2  +  (KAz(0x,0v)0x)2 


cos 


KA,tMy>e»l  \ 
m(0x,9y)  j  / 


+  WAz+<  1  > 


(8) 
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This  is  given  by  Az+  =  (Sum  +  Az)/Vl 
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For  small  tracking  errors  (i.e.,  191  <  X/d)  the  trigonometric  functions  can  be  linearized  and  Az  is 
linear  in  8  ."  In  order  to  circumvent  the  dependence  of  NEA  on  tracking  loop  bandwidth,  we  will  focus 
on  the  (single-sided)  noise  equivalent  spectral  density  (NESD).  The  quantities  are  related  by 


NESD  = 


NEA 

NEB 


(9) 


where  NEB  is  the  single-sided  noise  equivalent  bandwidth  of  the  closed-loop  tracking  system.  It  is  easy 
to  show  that  the  NESD  is  equal  to 


NESD  =  |~^K.  (0)21  1 
hv  Az 


(10) 


NESD  = 


4  (  ff  S(x,y)  L{x,y)  kx  dx  dy  )2 

0  0 


hv 


oooo 

JJ I  L(x,y)  I2  dx  dy 

0  0 


-1 


(11) 


From  the  Schwartz  inequality  the  optimal  local  oscillator  field  distribution  and  minimum  single-axis 
NESD  arc  equal  to 


L(x.y)opt  =  -S(x,y) 


(12) 


NESDop{  = 


-ooo© 


(13) 


For  tracking  in  elevation,  xS(x,y)  would  be  replaced  with  yS(x,y). 

As  might  be  expected  the  optimal  local  oscillator  field  distribution  is  the  product  of  the  distance 
from  the  center  of  the  detector  and  the  signal  field  amplitude.  The  distance  term  can  be  thought  of  as  a 
"lever  arm."  where  a  longer  lever  arm  implies  more  sensitivity  to  angular  information.  Note  that  for  an 
infinitely  large  and  infinitely  fine  detector  array  the  local  oscillator  could  be  arbitrarily  shaped  and  with 
the  proper  linear  processing  of  each  detector  element  output  the  same  bound  can  be  achieved.  To 
simultaneously  achieve  both  azimuth  and  elevation  bounds  would  require  such  processing,  since  one 
local  oscillator  field  cannot  simultaneously  optimally  process  azimuth  and  elevation.  It  is  easy  to  show 
that  for  two  dimensional  tracking,  an  infinite  area  quadrant  detector,  and  equal  weighting  of  azimuth  and 
elevation  tracking  errors  L(x,y)  =  (x+y)S(x,y)  best  minimize  the  overall  tracking  error.  Other  suboptimal 
local  oscillator  fields  will  be  discussed  later. 


“  Note  that  the  sign  of  the  discriminator  is  dependent  on  the  sign  of  Ao>.  If  there  is  uncertainty  in  Ato,  this 
must  be  taken  into  account  in  order  to  maintain  tracking-loop  stability. 
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It  was  previously  stated  that  the  optimal  tracking  performance  bounds  should  be  the  same  in  the 
focal  plane  as  the  pupil  plane.  We  will  show  that  the  respective  bounds  on  NESD  are  the  same.  It  is 
well  known  that  under  the  paraxial  propagation  approximation,  the  field  distribution  in  the  focal  plane 
(or  far  field)  can  be  related  to  the  field  distribution  in  the  pupil  plane  by  a  scaled  two  dimensional 
Fourier  transform  multiplied  by  a  phase  term  [7].  Using  this  property  as  welt  as  ParscvaTs  theorem  the 
above  expressions  can  be  rewritten  as 

L(X’y)  =  ^|“  (,4) 

NESDop.=  [i~  f  T~-y>  2  dx  “y]  '*  (,5) 


where  S(x,y)  and  L(x, y)  are  the  signal  and  local  oscillator  amplitudes  in  the  focal  plane,  respectively. 
This  bound  is  die  same  as  that  derived  for  the  noise  equivalent  spectral  density  of  a  heterodyne  spatial 
tracking  system  operating  in  the  focal-plane  [1].  Note  that  it  is  3  dB  worse  than  the  (minimum  mean- 
square  estimation)  bound  for  a  signal  shot-noise  limited  direct  detection  spatial  tracking  system[8].  It  is 
not  surprising  that  the  heterodyne  detection  performance  bounds  in  the  focal  plane  and  pupil  plane  are 
the  same  since  the  two  field  distributions  are  related  by  a  linear  Fourier  transform.  As  one  might  expect, 
the  two  optimal  processing  techniques  (LO  shapes)  are  also  related  by  a  Fourier  transform.  The  easy 
access  to  both  the  pupil  plane  and  its  transform  the  focal  plane  is  an  attractive  feature  of  optical  systems. 
However,  these  two  processing  techniques  are  analogous  to  amplitude  comparison  monopulse  and  phase 
comparison  monopulse  which  are  used  in  conventional  radar  systems  [9,10]. 

In  the  pupil  plane,  plane  wave  and  Gaussian  field  distributions  can  be  represented  by 


Plane  Wave 


Gaussian 


(16) 

(17) 


where  d  is  the  diameter  of  the  plane  wave,  o)$  is  the  1/e  amplitude  of  the  Gaussian  field,  and  circ(x)  is  1 
for  Ixl  <  1  and  0  otherwise.  The  focal-plane  amplitudes  are  then  given  by 


S(x,  y) 


2 J.  (M/717) 

ltd  1  Uf  ’  ' 


Xf  y2 


Airy  Disc 


(18) 
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where  f  is  the  focal  length  of  the  lens  used  to  transform  from  the  pupil  plane  to  the  focal  plane. 

To  include  the  effect  of  tracking  errors  one  either  multiplies  S(x,y)  by  exp[-jk(x0x  +  yO  )]  in  the 
pupil  plane  or  shifts  S(x,y)  by  S(x  -  fB^.y  -  f0y)  in  the  focal  plane.  The  resulting  bounds  on  tne  NESD 
are  given  by 

NESD  >.  f  ^ —  1  1  Plane  Wave  or  Airy  Disc  r 

I  hv  4  ^  2J 


NESD> 

lhv  &)2 


-i 


Gaussian 


(21) 
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3.  SUBOPTIMAL  LOCAL  OSCILLATORS  DISTRIBUTIONS 


The  on-axis  noise-induced  tracking  bounds  developed  in  the  previous  section,  which  are  valid  for 
an  infinitely  large  and  infinitely  fine  detector  array,  were  shown  to  be  independent  of  the  detection  plane. 
We  also  note  that  it  is  well  known  that  the  bound  on  the  sum  channel  mode-matching  efficiency  is  also 
independent  of  the  detection  plane.  In  this  section,  the  detector  array  will  be  restricted  to  be  a  simple 
quadrant  detector  and  the  local  oscillator  will  be  restricted  to  one  that  achieves  the  same  performance  in 
azimuth  and  elevation.  Although  under  such  restrictions  for  a  given  local  oscillator  pattern  the  sum 
channel  mode-matching  efficiency  is  still  independent  of  the  detection  plane,  the  tracking  bounds  are  not. 


We  will  investigate  two  suooptimal  local  oscillator  patterns.  A  plane-wave  signal  field  is  assumed 
since  it  closely  models  the  received  signal  for  most  systems.  Two  suboptimal  local  oscillators  to  be 
considered  are  a  plane  wave  matched  to  the  signal  field  and  a  truncated  Gaussian  field.  A  plane-wave 
local  oscillator  is  of  interest  since  it  will  maximize  the  power  in  the  sum  channel  thus  maximizing 
communication  or  radar  signature  performance.  A  truncated  Gaussian  field  is  of  interest  because  many 
lasers  naturally  emit  Gaussian  fields,  which  cannot  be  precisely  shaped  into  a  plane  wave  without  a 
prohibitive  loss  in  the  useable  local  oscillator  power.  For  a  plane-wave  local  oscillator  field  matched  to 
the  signal 


NESD  = 


l  hv  9  i^xyl 


-t 


Matched  Plane  Wave  in  Pupil  Plane 


(22) 


Thus  for  a  plane-wave  signal  field,  a  local  oscillator  mode  matched  to  the  sum  channel  is  within  1.4  dB 
of  the  optimal  single-axis  tracking  performance  (Equation  20).3  It  is  of  interest  to  compare  this  to  the 
focal-plane  case  where  the  local  oscillator  is  matched  to  the  Airy  disc.  This  was  shown  to  be  [  1  ] 

NESD  =  ( i  1  Matched  Airy  Disc  in  Focal  Plane 

'  9K  (?)  (2 

Thus,  in  the  focal  plane  a  local  oscillator  matched  to  the  signal  field  incurs  over  5.3  dB  of  loss  with 
respect  to  the  single-axis  tracking  bound.  The  goal  of  simultaneously  optimizing  communication  and 
tracking  performance  with  a  quadrant  detector  is  more  closely  met  in  the  pupil  plane  than  in  the  focal 
plane. 


3 In  this  paper,  a  comparison  in  terms  of  dB  refers  to  the  increase  in  signal  power  that  is  required  in  order 
to  maintain  the  same  NESD, 
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In  the  pupil  plane,  if  the  signal  is  a  plane  wave  and  the  local  oscillator  is  a  truncated  Gaussian  field, 
then  the  sum  channel  mode  matching  efficiency  and  the  pupil  plane  and  focal  plane  NESD  are  given  by 
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Pupil  Plane  (25) 
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Focal  Plane 
(26) 


where  o>j  is  the  1/e  amplitude  of  the  local  oscillator  field.  These  three  expressions  are  plotted  in  Figure  5 
relative  to  their  optimum  values.4  In  the  limit  as  (d/t»>j)  becomes  small  it  is  straightforward  to  show  that 
these  results  agree  with  the  matched  plane  wave  and  Airy  disc  results  above.  For  d/tOj  <  2.3  the  tracking 
performance  of  the  pupil-plane  system  is  superior.  Ford/Wj  >  2.3  the  focal-plane  system  is  superior. 

This  is  expected  since,  in  the  pupil  plane,  the  most  sensitive  angular  information  is  near  the  edges  of  the 
signal  field.  As  the  local  oscillator  field  becomes  more  and  more  Gaussian,  this  angular  information  is 
"amplified"  by  the  local  oscillator  field  less  and  less.  Conversely,  in  the  focal  plane  as  the  local  oscilla¬ 
tor  field  becomes  more  Gaussian,  its  diffraction  pattern  in  the  focal  plane  becomes  wider.  At  first  this 
wider  local  oscillator  pattern  in  the  focal  plane  increasingly  weights  the  angular  sensitive  area  of  the 
signal  pattern  (areas  of  large  slope)  until  an  optimum  value  of  d/o>j  =  3.1  is  reached.  Beyond  this  value 
the  local  oscillator  pattern  is  spread  too  widely  and  the  performance  suffers. 

In  order  to  maximize  the  mode  matching  efficiency  in  the  sum  channel,  the  desire  to  shape  the  local 
oscillator  pattern  into  a  plane  wave  is  clear.  Of  course  the  power  in  the  local  oscillator  field  delivered  to 
the  detectors  is  reduced  as  the  beam  is  shaped.  This  power  loss  is  simple  to  calculate  and  is  given  by 

LO  Power  loss  =  1  -  exp(  -  ~  (^-)2 )  .  (27) 

This  expression  is  also  plotted  in  Figure  5.  There  is  a  large  penalty  in  delivered  local  oscillator  power  for 

d/to)j  <  1 . 


4 


The  optimum  value  for  m  is  1  and  the  optimum  value  for  both  NESD  values  is  given  in  Equation  20. 
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RATIO  OF  SIGNAL  DIAMETER  TO  LO  1/a  RADIUS  (  £  ) 

“l 

Figure  5.  Comparison  of  pupil  plane  and  focal  plane  sum  and  difference  channel  mode-matching  efficiencies  and 

LO  power  loss. 


As  an  example  of  the  trade-off  between  simultaneously  optimizing  communication  and  tracking 
performance  we  refer  to  the  Laser  Intersatellite  Transmission  Experiment  (LITE).  The  anticipated 
operating  point  for  the  LITE  local  oscillator  field  is  d/w(  =  2  [Reference  11].  This  value  was  chosen  to 
minimize  the  communication  channel  mode  matching  loss  (<  0.35  dB)  while  also  minimizing  the  loss  in 
local  oscillator  power  (<  0.63  dB).  The  desire  to  minimize  the  communication  signal  power  loss  is 
obvious.  The  need  to  minimize  the  loss  of  local  oscillator  power  stems  from  the  desire  to  keep  the  signal 
levels  above  the  thermal  noise  in  the  front-end  amplifiers.  At  this  value  of  d/u>t  pupil-plane  processing 
only  slightly  out-performs  the  focal-plane  system  (~  0.8  dB).  The  pupil-plane  tracking  performance  is 
approximately  2.7  dB  from  optimum  and  the  focal-plane  tracking  performance  is  approximately  3.5  dB 
from  optimum. 

So  far  we  have  only  considered  coherent  detection  followed  by  coherent  demodulation  of  the  I F 
currents.  In  practice,  due  to  nonzero  laser  line-widths  and  low  received  signal  power,  coherent  demodu¬ 
lation  is  not  possible.  For  high  I.F.  SNR  ratios  the  NESD  expressions  are  still  valid.  However,  for  low 
I F  SNR  the  expressions  for  tracking  performance  arc  worse  than  these  equations  predict  [1].  The  exact 
details  of  the  non-coherent  demodulation  penalty  is  very  dependent  on  the  particular  processing.  How¬ 
ever,  for  many  demodulation  schemes  (i.e.,  square-law  and  envelope)  the  NESD  will  be  proportional  to 
the  product  of  the  difference  channel  mode  matching  efficiency  and  the  sum  channel  mode  matching 
efficiency.  That  is 
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The  local  oscillator  distribution  (for  two-axis  tracking  and  assuming  a  quadrant  detector)  which  maxi¬ 
mizes  these  expressions  is  easily  found  by  optimizing  the  linear  combination  of  either  { S(x,y),  dS(x, y)/ 
dx+  5S(x,y)/dy )  for  focal-plane  processing  or  { S(x,y),  xS(x,y)  +  yS(x,y) }  for  pupil-plane  processing.  It 
was  previously  found  that  in  the  case  of  focal-plane  processing,  that  (0.3  IS)2  bounded  the  quantity  in 
square  brackets5.  For  pupil-plane  processing  the  quantity  in  square  brackets  is  less  than  (0.40 1)2. 

Thus,  there  is  only  1.0  dB  of  improvement  available  with  pupil-plane  tracking  performance  when  operat¬ 
ing  on-axis  at  low  SNR.  It  was  also  determined  that  a  local  oscillator  distribution  matched  to  the  sum 
channel  ((0.2 12)2)  was  1.7  dB  from  its  low  SNR  bound  [1].  For  pupil-plane  processing  a  local  oscillator 
distribution  matched  to  the  sum  channel  [(0.333)2]  is  within  0.80  dB  from  this  bound. 

The  performance  of  the  pupil-plane  and  focal-plane  tracking  systems  operating  at  low  I F  SNR 
with  a  local  oscillator  that  is  a  truncated  Gaussian  field  in  the  pupil  plane  can  be  seen  in  Figure  5.  The 
relative  tracking  performance  (in  dB)  is  given  by  one-half  the  distance  between  the  pupil-plane  and  the 
focal-plane  tracking  curves.  For  instance,  for  d/tOj  =  0.0  the  pupil-plane  tracking  system  is  superior  by 
2.0  dB  (this  is  consistent  with  the  previous  paragraph),  for  d/u>j  =  2.3  the  pupil-plane  and  focal-plane 
systems  have  the  same  performance,  and  for  d/tu.  =  4.0  the  focal-plane  tracking  system  is  superior  by 
2.7  dB.  1 


5In  Reference  [1]  the  quantity  is  listed  as  0.2832.  The  difference  stems  from  the  fact  that  in  [1]  the  radial 
derivative  was  used  instead  of  dS/dx  +  dS/dy. 
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4.  OFF-AXIS  DISCRIMINATOR  GAIN 


As  we  have  seen,  a  quadrant  detector  can  come  very  close  to  achieving  optimal  tracking  perform¬ 
ance.  We  now  investigate  some  of  the  issues  associated  with  off-axis  discriminator  gain  (a  major  factor 
affecting  pull-in  performance)  of  a  quadrant  detector  based  tracking  system.  It  is  straightforward  to 
show  that  the  sum  channel  discriminators  in  the  focal  olanc  and  the  pupil  plane  are  identical  by  invoking 
the  Fourier  transform  relationship.  However,  this  is  not  necessarily  true  for  the  difference  channels. 
Although  the  slope  around  null  of  the  difference  channel  angle  discriminators  may  be  the  same  in  the 
focal  plane  and  the  pupil  plane,  their  gain  away  from  null  is  not  necessarily  the  same.  This  is  not  to  say 
that  no  local  oscillator  distribution,  electrical  processing,  and  detector  will  yield  precisely  the  same 
discriminator  profiles  in  the  focal  plane  and  pupil  plane,  but  rather  that  for  a  quadrant  detector  the 
discriminators  are  not  necessarily  the  same.  To  illustrate  this  fact,  the  expressions  for  pupil-plane  and 
focal-plane  discriminators  with  Gaussian  signal  and  local  oscillator  distributions  are  given  by 
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These  are  plotted  as  a  function  of  full-width  half-power  bcamwidths  in  Figure  6  (one  beamwidth  =  2.355/ 
(k<i)s)  =  0.375X/o)s).  Around  boresight  the  discriminators  have  the  same  slope  and  hence  the  same  track¬ 
ing  performance  (this  is  in  contrast  to  the  previous  plane-wave  signal  and  local  oscillator  fields  which 
yielded  discriminator  gains  that  were  different  by  dB).  As  0  becomes  large,  the  focal-plane  discrimina¬ 
tor  falls  off  as  exp[-(k0a>s/V8)2],  whereas  the  pupil-plane  discriminator  falls  off  only  as  l/(k0(os/V8).  For 
this  case,  the  pupil-plane  tracking  system  has  the  potential  of  having  much  better  pull-in  performance. 

We  now  investigate  the  focal-plane  and  pupil-plane  discriminator  properties  for  a  plane-wave 
signal  pattern  and  a  Gaussian  local  oscillator  pattern  which  has  been  truncated  in  the  pupil  plane.  Equa¬ 
tions  for  the  angle  discriminator  profiles  are  given  below 
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ANGULAR  OFFSET  (BMmwkfths  kOu>/2.355) 

Figure  6.  Comparison  of  matched  Gaussian  focal-plane  and  pupil-plane  discriminators. 
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where  sgn(x)  is  -1  for  x  <  0  and  1  for  x  >  0.  These  equations  are  plotted  as  a  function  of  the  diffraction 
limited  beamwidth  (X/d)  for  five  different  ratios  of  d/w.  ( 0. 0.2,  2. 4,  20)  in  Figures  7  to  1 1 .  The  corre¬ 
sponding  on-axis  sum  channel  power  losses  are  (0, O.Oi,  0.34. 3.17, 17  dB).  respectively.  The^c  Figures 
are  consistent  with  the  on-axis  sum  and  difference  gams  shown  in  Figure  5.  The  case  where  d/o>(  =  0 
corresponds  to  matched  signal  and  local  oscillator  patterns.  This  represents  the  optimal  sum  channel 
mode-matching  efficiency  and  as  was  stated  earlier,  the  gain  in  the  linear  region  is  ~  4  dB  larger  in  the 
pupil  plane  than  the  focal  plane.  Note  that  the  gam  of  the  difference  channel  away  from  null  in  the 
pupil  plane  only  slightly  out-performs  the  gain  in  the  focal  plane.  There  is  also  a  large  amount  of 
ringing  in  the  difference  and  sum  channels  as  compared  with  the  Gaussian  case.  The  local  oscillator 
distribution  can  be  shaped  to  decrease  this  ringing  and  increase  the  off-axis  gain  at  the  expense  of  gain  in 
the  linear  region  and  sum  channel  mode  matching  efficiency  as  is  seen  in  the  figures.  A  strong  analogy 
exists  between  this  effect  and  the  trade-off  of  side-lobe  suppression  and  on-axis  gain  that  is  well  known 
in  conventional  antenna  design  (12]. 


These  figures  reinforce  the  earlier  statement  that,  at  both  high  and  low  I F  SNR.  for  d/u»(  <  2.3  the 
pupil-plane  system  has  superior  discriminator  gain,  while  for  d/u>j  >  2.3  the  focal-plane  is  superior. 
However,  note  that  in  almost  ail  cases  for  large  off  axis  angles  the  pupil-plane  system  has  superior 
difference  channel  discriminator  gain.  Therefore,  the  pupil-plane  systems  have  the  potential  for  superior 
pull-in  performance. 
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Figure  #.  Comparison  of  focal-plane  and  pupil-  plane  discriminators. 
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Figure  9.  Comparison  of  focal-plane  and  pupil-plane  discriminators. 
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5.  APPLICATIONS  OF  PUPIL-PLANE  PROCESSING 
USING  HIGHER  ORDER  ARRAYS 


We  have  mainly  focused  on  discussing  results  applicable  to  a  spatial  tracking  system  which  employ 
a  quadrant  detector  in  the  pupil  plane.  In  this  section,  we  briefly  comment  on  the  potential  applications 
of  using  higher  order  detector  arrays.  As  was  stated  previously,  almost  all  of  these  are  applicable  to  the 
focal  plane;  however,  implementation  considerations  seem  to  favor  the  pupil-plane  approach.  Note  that 
almost  all  the  applications  discussed  have  direct  analogies  in  phased-array  radars. 

The  use  of  an  N  x  N  detector  array  allows  the  receiver  FOV  to  be  electrically  scanned  by  adding  a 
phase  gradient  across  the  detector  array  outputs.  This  makes  possible  high  speed  scanning  for  acquisi¬ 
tion  as  well  as  increased  pull-in  performance.  The  limit  on  the  FOV  is  when  the  tilt  across  a  detector 
element  approaches  its  diffraction  limited  beamwidth,  i.e.. 


0 


max 


(38) 


Beyond  this  point,  the  mode  matching  efficiency  within  each  detector  element  decreases.  Thus,  for  a 
10X10  detector  array  the  FOV  can  be  scanned  or  the  pull-in  increased  to  -10  beamwidths. 

Such  a  detector  array  could  also  provide  a  point-ahead  function  for  a  lasercom  link,  thus  eliminat¬ 
ing  the  bulky  point-ahead  device.  By  boresighting  the  tracking  detector  with  no  phase  gradient  and  then 
adding  a  phase  gradient,  a  look-behind  angle  can  be  incorporated  into  the  tracking  system.  The  point- 
ahead  requirement  for  a  GEO  to  LEO  link  is  -  40  p,rad,  Assuming  a  typical  antenna  beamwidth  of 
4  firad,  a  10  x  10  detector  array  would  be  required  to  fulfill  the  point-ahead  requirement. 

Atmospheric  compensation  could  also  be  performed  using  pupil-plane  processing  [13],  It  is  well 
known  that  when  communicating  through  the  atmosphere,  if  the  ratio  of  antenna  diameter  to  transverse 
coherence  length  (d/rQ)  is  much  larger  than  1,  then  there  is  a  stiff  mode  matching  penalty  (>10  dB).  This 
penalty  could  be  eliminated  by  servoing  the  phases  of  adjacent  detectors  in  multiple  detector  array  placed 
in  the  pupil  plane  provided  that 


Note  that  the  well-known  alternative  of  using  a  deformable  mirror  has  the  advantage  of  also  compensat¬ 
ing  the  transmitted  beam  as  well  as  the  received  beam. 

Other  possible  applications  include  multiple  receiver  apertures,  adaptive  antenna  nulling,  and  lateral 
and  longitudinal  mode  stabilization  of  laser  diode  arrays. 


21 


6.  EXPERIMENTAL  RESULTS 


The  optical  layout  and  signal  processing  electronics  are  illustrated  in  Figures  12  and  13.  Most  of 
the  experimental  setup  and  electronics  were  originally  designed  for  use  in  a  focal-plane  tracking  sys¬ 
tem  [141.  Consequently,  the  optics  and  electronics  have  not  been  completely  optimized  for  a  pupil-plane 
tracking  system.  However,  this  system  demonstrates  a  number  of  key  features  of  a  single-axis,  closed- 
loop.  heterodyne  spatial  tracker  using  pupil-plane  processing. 

The  tracker  was  organized  as  an  integral  part  of  a  balanced  communications  receiver.  Separate 
tracking  and  communication  detectors  are  used.  Throughout  the  discussion  of  the  single-axis  experi¬ 
ment,  track  channel  references  apply  only  to  the  azimuth  channel.  The  track  channel  power,  P  ,  is 
defined  as  the  signal  optical  power  in  the  azimuth  track  channel  before  the  image  splitter,  or,  equival¬ 
ently,  as  the  sum  of  the  power  at  the  right  and  left  track  detectors.  The  communications  channel  power, 

P  .  is  defined  as  the  sum  of  the  signal  optical  power  at  both  communications  detectors.  The  input  signal 
power,  P  ,  refers  to  the  total  signal  optical  power  at  the  input  to  the  combining  beamsplitter  (for  lossless 
optics),  and  is  equal  to  P  +  P  . 

The  signal  beam  was  passed  through  galvonometer-driven  disturbance  and  tracking  mirrors,  and 
then  combined  with  the  local  oscillator  beam  on  a  50:50  (s  and  p  polarization)  combining  beamsplitter. 
The  signal  and  local  oscillator  lasers  were  GaAlAs  semiconductor  diode  lasers  operating  at  a  wavelength 
of  0.86  p.m.  The  optical  beams  were  each  collimated  Gaussian  beams  with  1/e  radius  =  2.65  mm. 

The  rms  wavefront  quality  was  interferometrically  measured  to  be  A/33  and  X/29  (or  equivalently  Strehl 
ratios  of  0.92  and  0.90)  for  the  local  oscillator  and  signal  beams,  respectively.  The  overlapped  signal  and 
local  oscillator  beams  were  directed  to  the  intensity-noise-cancelling  communications  detectors.  The 
sum  channel  mode  matching  efficiency  was  measured  to  be  0.85  using  the  technique  described  in  [15]. 
The  polarization  beatr.  plitters  in  each  of  the  communications  paths  were  used  to  direct  a  portion  of  the 
received  signal  power  to  the  tracking  optics.  In  the  single-axis  demonstration  only  the  azimuth  tracking 
channel  was  actually  used.  The  tracking  beam  was  split  by  a  knife  edge,  from  which  the  light  was 
focussed  onto  two  tracking  detectors. 

The  polarization  of  the  signal  and  local  oscillator  beams  at  the  input  to  the  combining  beamsplitter 
was  separately  adjusted  (using  the  half-wave  plates)  such  that  equal  local  oscillator  power  was  directed 
to  all  four  detectors,  while  the  signal  power  ratio,  P  /P  ,  was  set  to  81/19  between  the  communications 
and  track  paths.  The  signal  power  split  represented  a  compromise  between  providing  as  much  communi¬ 
cations  pov/er  as  possible  and  still  maintaining  a  sufficiently  robust  track  signal  such  that  adequate  track 
performance  was  achieved. 

Just  ahead  of  the  combining  beamsplitter  a  quarter-wave  plate  was  placed  in  the  signal  path,  with  its 
axes  aligned  with  the  polarization  beamsplitters'  axes.  This  introduced  a  90°  phase  shift  between  the 
communication  I  F  signal  and  each  of  the  track  channel  I  F  signals.  Such  a  phase  shift  rendered  the  sign 
of  the  demodulated  error  signal  insensitive  to  changes  in  the  sign  of  Au>  (see  Footnote  2),  the  difference 
between  the  local  oscillator  und  signal  laser  frequencies.  In  other  words,  the  sign  of  the  resultant  dis¬ 
criminant  was  not  affected  by  whether  the  local  oscillator  operates  above  or  below  the  signal  laser  in 
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Figure  12.  Optical  layout,  I -axis. 
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Figure  13.  Correlation  processing. 


frequency.  We  note  that  the  phase  shift  must  be  introduced  optically  and  not  electrically  in  order  to 
accomplish  this  effect.  The  need  to  render  the  tracking  error  signal  insensitive  to  Ao>  may  or  may  not  be 
necessary  in  a  particular  system.  It  is  dependent  on  frequency  uncertainty  during  acquisition  and  the 
frequency  stability  during  tracking. 


The  outputs  from  the  tracking  and  communications  detectors  were  amplified  and  combined  as 
shown  in  Figure  13.  Ideally,  one  would  use  a  truly  coherent  processing  scheme,  which  would  yield  an 
NEA  given  by  (Equation  10) 


NEA 


ideal 


ka!<°> 


(40) 


where  y  is  included  to  account  for  the  nonideal  phase  and  amplitude  distortions  that  give  rise  to  a  Strehi 
ratio  less  than  1.  For  matched  Gaussian  beams  KAz(0)  =  0.94  BW'1.  Unfortunately  coherent  processing 
is  quite  difficult  to  implement  unless  either  the  signal  frequency  is  well  known  or  the  I  F  SNR  is  much 
greater  than  one.  Because  the  heterodyne  tracker  is  designed  to  work  in  a  communications  system  which 
employs  FSK  modulation  over  a  500  MHz  bandwidth  with  low  received  signal  powers,  neither  of  these 
conditions  is  met.  Two  alternatives  include  either  square-law  power  detection  of  the  heterodyne  power 
in  each  side  of  the  track  channel,  or  the  correlation  technique  illustrated  in  Figure  13. 

Square-law  detection  suffers  from  excessive  sensitivity  to  noise-induced  bias  at  low  input  SNR, 
where  the  dc  term  due  to  squared  broad  band  noise  may  be  10  to  50  times  the  dc  term  due  to  the  signal 
power.  This  dc  noise  term  is  successfully  removed  during  subtraction  of  the  right  and  left  I  F  channel 
signal  powers  only  if  the  right  and  left  channels  are  phase  and  amplitude  matched  to  better  than  -  5  % 
over  the  full  input  bandwidth,  which  is  difficult. 

An  alternative  correlation  scheme  employed  in  this  experiment  used  the  I  F  communications  signal 
as  a  coherent  reference  to  shift  the  track  signals  to  baseband  without  producing  a  dc  bias,  as  long  as  the 
communications  and  track  channel  noises  were  uncorrelated  [  14],  Prior  to  correlation  the  communica¬ 
tions  and  track  channels  were  offset  in  frequency  (via  SSB  upconversion)  by  signals  which  differ  in 
frequency  by  83  kHz.  The  resultant  post-correlation  signals  were  centered  at  this  difference  frequency, 
instead  of  at  baseband,  in  order  to  reduce  sensitivity  to  the  dc  offsets  of  the  wideband  mixers.  Subse¬ 
quent  synchronous  demodulation  shifted  the  error  signals  to  baseband.  Baseband  processing  produced 
an  error  signal  normalized  by  the  total  signal  power.  Normalization  was  provided  in  order  to  maintain 
constant  error  gain,  and  hence  constant  closed-loop  bandwidth,  over  a  wide  range  of  optical  signal 
powers. 

The  correlation  scheme  used  to  generate  the  error  discriminant  is  sensitive  to  any  noise  which  is 
common  to  the  communications  and  track  channels.  Such  noise,  for  instance,  is  caused  by  intensity 
fluctuations  in  the  local  oscillator  laser.  Unlike  the  shot  and  amplifier  thermal  noises,  which  are  uncorre¬ 
lated  between  ali  of  the  detectors,  uncancelled  intensity  noise  can  cause  an  offset  error  in  the  tracker 
which  may  limit  performance.  However,  the  balanced  receiver  1 16]  used  in  the  communications  channel 
achieved  sufficient  cancelling  (greater  than  20  dB)  to  make  the  bias  due  to  uncancelled  intensity  noise  a 
small  fraction  of  a  beamwidth. 
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For  a  shot-noise-limited  system  the  idea!  NEA  for  the  correlation  scheme  can  be  written  as 


NEA 


_  1 


NEB  r 


i<fea'  KAz<0)  x 

hv 


I  + 


_w_- 
mfc 
hv  - 


1/2 


(41) 


where  W  is  the  noise  equivalent  bandwidth  of  the  i.f.  bandpass  filters  The  first  two  terms  are  identical  to 
the  ideal  coherent  demodulated  NEA  result  given  in  Equation  40.  The  term  in  the  square  brackets  can  be 
interpreted  as  a  noncoherent  detection  penalty  and  is  almost  identical  to  that  given  in  [ !  ].  Note  that  for 
high  1 F  SNR  there  is  no  penalty  for  noncoherent  detection. 

The  measured  experimental  results  were  worse  that  the  ideal  result  given  in  Equation  41  for  several 
reasons.  First,  excessive  thermal  noise  in  the  detector  amplifiers  made  it  difficult  to  supply  enough  local 
oscillator  power  to  the  detectors  to  achieve  the  snot-noise  limit.  The  communications  channel  noise  was  a 
factor  of  2  above  shot  noise,  while  the  track  channel  operated  a  factor  of  3  above  shot  noise.  Secondly,  as 
noted  above,  the  particular  electronic  processing  used  (SSB  upconversion)  in  .'his  experiment  doubled  the 
expected  noise  level  (\2  increase  in  NEA).  The  resultant  expected  experimental  NEA  can  be  written  as 


NEA 


experiment 


Jl 

Ka/°> 


(42) 


The  unnormalized  discriminant  was  measured,  and  is  compared  with  the  predicted  discriminant  in 
Figure  14.  The  predicted  discriminant  for  the  correlation  based  tracker  is  the  product  of  the  sum  channel 
and  difference  channel  discriminator  profiles  shown  in  Figure  6.  The  on-axis  gain  was  within  2%  of  the 
theoretically  predicted  gain  for  matched  Gaussian  beams  (Equations  30  and  31),  while  the  off-axis  behav¬ 
ior  is  also  in  excellent  agreement  with  theory.  The  communications  channel  falls  off  in  amplitude  more 
quickly  than  either  half  of  the  track  channel,  which  means  that  the  off-axis  performance  of  the  communi- 
cations-track  product  is  dominated  by  the  communications  channel. 

The  rms  noise  over  a  bandwidth  of  1  to  3  kHz  was  measured  at  the  tracker  output  over  a  25  dB 
range  of  input  signal  optical  power.  In  order  to  specifically  characterize  the  shot-noise  limit  of  the  angle 
sensor,  the  noise  was  measured  with  the  feedback  loop  open,  and  at  frequencies  well  above  the  acoustic 
and  mechanical  disturbances  which  dominated  the  noise  spectrum  below  200  Hz  with  the  loop  open. 
Because  the  shot  noise  spectrum  was  flat  within  the  measurement  region  of  interest  (0  to  3  kHz),  the  shot 
noise  contribution  within  a  2  kHz  bandwidth  was  the  same  from  1  to  3  kHz  as  from  0  to  2  kHz.  The  NEA 
was  obtained  by  dividing  this  rms  noise  by  the  measured  discriminant  gain  at  boresight.  The  resultant 
NEA  vs  input  optical  power  for  the  pupil-plane  tracker  is  plotted  in  Figure  15.  In  addition,  the  NEA 
bound  using  the  optimum  local  oscillator  (Equation  21)  and  tj  -  1,  the  ideal  NEA  produced  via  correla¬ 
tion,  with  a  matched  Gaussian  local  oscillator,  -q  =  0.73,  and  y  =  0.85,  and  the  predicted  experimental 
NEA  are  plotted  for  comparison.  We  note  that  the  measured  NEA  was  1.7  dB  (in  required  input  signal 
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MEASURED  DISCRIMINANT  OUTPUT  FOR 
FIXED  OFFSET  Of  ♦  0.23  BW 


MEASURED  DISCRIMINANT  OUTPUT  FOR 
FIXED  OFFSET  OF  -  0.23  BW 
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figure  14.  Pupibplane  discriminant. 
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Figure  15.  NEA  vs  signal  optical  power. 
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power)  above  the  predicted  (experimental  model)  NEA.  Further  measurements  isolated  this  discrepancy 
to  the  communications  channel  SNR  prior  to  correlation.  If  the  measured  SNR  in  the  communication 
channel  is  used  to  predict  the  experimental  NEA,  excellent  agreement  results,  which  validates  the 
analysis  of  the  correlation  process. 

The  -3  dB  closed  loop  bandwidth  of  the  tracking  loop  was  measured  at  1.3  kHz.  Higher  bandwidth 
could  have  been  achieved  but  was  not  necessary  for  this  application  [81.  The  closed  loop  and  rejection 
transfer  functions  are  shown  in  Figure  16.  Normalization  of  the  on-axis  gain  by  the  sum  channel  signal 
made  the  closed  response  of  the  tracker  insensitive  to  input  signal  optical  power  variation  over  a  23  dB 
range. 
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(a)  CLOSED  LOOP  TRANSFER  FUNCTION 


10  Hz  100  Hz  1kHz 
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(b)  REJECTION  TRANSFER  FUNCTION 
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Figure  16  Heterodyne  tracker  closed  loop  response. 
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7.  CONCLUSION 


We  have  shown  that  the  theoretical  tracking  bounds  for  a  spatial  tracking  system  using  coherent 
detection  in  the  pupil  plane  is  identical  to  the  more  conventional  focal-plane  processing  technique.  More 
important,  we  have  shown  that  when  implementation  issues  such  as  receiver  complexity  are  considered, 
the  pupil-plane  system  offers  a  number  of  important  advantages  over  the  focal-plane  system.  For  in¬ 
stance.  with  a  quadrant  detector  the  optimum  communication  and  tracking  local  oscillator  fields  are  more 
closely  matched  in  the  pupil  plane  than  in  the  focal  plane.  This  fact  leads  to  a  potential  4  dB  improve¬ 
ment  over  focal-plane  processing.  Also  shown  was  an  increased  off-axis  discriminator  gain  which  leads 
to  increased  pull-in  performance.  Other  implementation  issues  such  as  detector  dead  zone,  cross  talk, 
and  alignment  sensitivity  were  shown  to  be  improved  in  the  pupil-plane  system. 

The  performance  of  a  Gaussian  local  oscillator  distribution  truncated  in  the  pupil  plane  was  also 
investigated.  It  was  shown  that  for  d/a>(  <  2.3  the  pupil-plane  tracking  system  was  superior  for  both  high 
and  low  I  F  SNR.  For  dA^  >  2.3  the  focal-plane  tracking  system  had  higher  discriminator  gain.  How¬ 
ever,  the  utility  of  this  region  is  minimized  by  the  fact  that  the  sum  channel  mode  matching  efficiency  is 
lower  making  this  region  an  undesirable  operating  point. 

Experimental  measurements  on  the  noise  performance  and  dynamics  of  a  single-axis  closed-loop 
tracking  system  using  detection  in  the  pupil  plane  were  presented.  Although  the  system  (especially  the 
electronics)  was  not  yet  optimized  for  pupil  plane  detection,  the  experimental  results  confirm  the  basic 
concept  and  many  of  the  theoretical  predictions. 
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